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STABILITY  PROPERTIES  OF  AZIMUTHALLY  SYMMETRIC 
PERTURBATIONS  IN  AN  INTENSE  ELECTRON  BEAM 

I.  INTRODUCTION 

In  recent  years  there  have  been  numerous  theoretical  investigations 
1  2  3-8 

of  the  equilibrium  '  and  stability  properties  of  intense  charged 

particle  beams,  motivated  by  a  variety  of  applications,  including  confine- 
9-11 

ment  and  transport  of  a  nonneutral  electron  beam,  inertial  confinement 
12  13 

fusion  '  driven  by  charged  particle  beams,  and  electron  beam  propaga- 
14  15 

tion  '  through  a  background  plasma.  Although  these  research  activities 
have  different  goals  and  objectives,  they  have  in  common  the  need  to  under¬ 
stand  the  stability  properties  of  intense  charged  beams  characterized  by 
strong  self  electric  and  magnetic  fields.  However,  previous  theoretical 
studies  of  the  resistive  instabilities  (except  for  hose)  have  been  based 
on  highly  simplified  models,  such  as  the  cold  fluid  treatment  of  Ref.  8. 

The  primary  purpose  of  the  present  paper  is  to  develop  a  Vlasov-Maxwell 
description  of  azimuthally  symmetric  perturbations  about  an  electron  beam, 
including  the  higher  radial  mode  numbers  (e.g.  hollowing  modes),  which 
have  been  given  scant  attention  in  the  literature,  although  they  appear  to 
be  troublesome  in  the  presence  of  a  return  current.  We  also  include  a 
treatment  of  the  hose  instability  based  on  the  same  model.  With  straight¬ 
forward  modifications,  the  results  of  this  paper  can  also  be  applied  to 
the  transverse  instabilities12'16,17  of  an  intense  ion  beam,  which  are 
particularly  important  in  heavy-ion  fusion  applications. 

The  equilibrium  and  stability  analysis  is  carried  out  for  an 
infinitely  long  electron  beam  propagating  parallel  to  a  uniform  applied 
magnetic  field  B^§  through  a  background  plasma  with  conductivity  o(r), 
which  in  general  can  be  a  function  of  the  oscillation  frequencies  of 

fluctuating  fields  and  the  collision  frequencies  of  each  species  in  the 
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plasma-  Our  solutions  are  carried  out  for  step-function  profiles  and 

real  values  of  a(r),  but  the  value  of  a  can  be  arbitrarily  large  or  small. 

The  radial  profiles  of  plasma  charge  and  current  are  assumed  to  be 

similar  to  the  beam  profile,  but  the  degree  of  charge  and  current 

neutralization  is  arbitrary.  It  is  assumed  that  v/Y^  <<  1,  where  v  and 
2 

Ybmc  are  Budker’s  parameter  and  the  characteristic  energy,  respectively, 
m  is  the  electron  rest  mass  and  c  is  the  speed  of  light  in  vacuo. 

Equilibrium  and  stability  properties  are  calculated  for  an  electron  beam 
described  by  a  "loss-cone"  distribution  function  [Eq.  (2) ],  which  leads 
to  a  flat-topped  beam  radial  profile.  Equilibrium  properties  and  basic 
assumptions  are  briefly  discussed  in  Sec.  II. 

The  formal  stability  analysis  for  azimuthally  symmetric  perturbations 
is  carried  out  in  Sec.  Ill,  within  the  framework  of  the  linearized  Vlasov- 
Maxwell  equations.  In  the  analysis,  the  response  of  the  background  plasma 
is  incorporated  in  terms  of  the  plasma  conductivity.  An  integro-differential 
eigenvalue  equation  [Eqs.  (22) ]  is  obtained  for  azimuthally  symmetric 
perturbations,  given  any  value  of  plasma  conductivity,  assuming  long  wave¬ 
length  (kR^  <<  1)  ,  low  frequency  (WR^  <<  c)  perturbations,  where  k  is  the 
axial  wavenumber,  uu  is  the  oscillation  frequency,  and  is  the  beam  radius. 
Equations  (22)  constitute  a  principal  result  of  this  paper  and  can  be  used 
to  investigate  stability  properties  for  a  broad  range  of  system  parameters. 

In  Sec.  IV,  dispersion  relations  for  the  radial  mode  numbers  n  *  0, 

1  and  2  are  obtained  analytically  from  the  integro-differential  eigenvalue 
equation  (22)  by  an  approximation  method  based  on  the  assumption  uut^  <  1, 


*  2  2 

where  td  is  a  magnetic  decay  time  defined  by  =  rc^R^^c  [Eq.  (26)  ]. 
For  example,  the  dispersion  relation  for  the  n  =  1  (sausage)  mode  can  be 
expressed  as  [Eq.  (65) j, 


z2  -  n  +  C*b  +  1Ab 


(5  -  i)  [l  -  i^UR^/c)  /6 J 


=  0, 


where  Z  *  (o»  -  is  the  normalized  Doppler-shifted  eigenfrequency, 

2  2  A  2  A 

n  =  2[0  (1  -  f)  -  (1  -  f  )  J  +  u>  , /a)  ,  and  £  =  4Tra./w  are  the  focussing 
Dm  e  cd  p d  1 

2  2  2 

force  and  conductivity  parameters,  respectively,  B,  *  (y.  -  l)/y.  , 

b  D  D 

A  2"  1/2 

u>  .  =  eB_/y_mc  and  &  ,  =  (47re  n,/y,  m)  are  the  beam  electron  cyclotron 
cd  u  b  pb  b  b 

and  plasma  frequencies,  respectively,  f  and  f  are  the  fractional  charge 

e  m 

and  current  neutralization  ,  respectively.  The  limiting  case  of  a  non¬ 
neutral  electron  beam  in  a  non-conducting  environment  (£-•  0)  is  discussed 
in  Sec.  IV.  For  example,  in  the  limit  of  5  *  0,  the  dispersion  relation 
for  n  =  1  [in  Eq.  (65) J  reduces  to 


z2  *  n  +  1/Y2  , 

which  is  identical  to  the  previous  result  obtained  by  Uhm  and  Davidson.3 

Analytic  investigations  of  the  dispersion  relations  for  n  =  0, 

1  and  2  are  carried  out  in  Sec.  V.A  for  collision-dominated  plasmas  with 

a  high  conductivity  satisfying  c,  »  1.  We  find  that  both  the  n  =  1 

(sausage)  and  n  =  2  (hollowing)  modes  are  driven  unstable  by  either  a 

strong  fractional  plasma  return  current  (f  )  or  by  a  magnetic  phase  lag 

m 

(ujt  >  0)  .  As  an  example  a  simplified  version  of  the  n  ■  1  dispersion 

a 

relation  [Eq.  (67) J, 


s'tu  - 


If  )  - 
m 


iC 

6 


CwR^/c)  ] 


3 


2 

which  is  obtained  from  Eq.  (65)  in  the  limit  of  n  =  26.  (1  -  f  )  and 

b  m 

<<  1,  shows  clearly  that  the  system  is  unstable  whenever  the 

fractional  current  neutralization  f  >0.5,  even  for  a  very  small  magnetic 

m 

decay  time.  It  is  also  remarkable  to  observe  that  for  f  =  0,  Eq.  (67) 

m 

8 

is  similar  to  the  result  obtained  by  Weinberg  within  the  framework  of  the 
coaxial  circular  orbit  beam  model. 

Several  points  are  noteworthy  in  the  numerical  analysis  of  the 
dispersion  relations  presented  in  Sec.  V.B.  for  a  self-pinched  electron 
beam  (u>  .  =  0) .  First,  the  n  =  1  perturbation  is  the  most  unstable 

CD 

azimuthally  symmetric  mode.  Even  if  f  =0,  the  maximum  growth  rate  of 

nt 

instability  for  the  n  =  1  perturbation  is  =  ImZ  =  0.205,  which  is 
comparable  to  that  of  the  resistive  hose  instability.7  Second,  the  n  =  2 
perturbation  has  two  unstable  modes.  Third,  except  the  n  *  0  perturbation, 
we  find  for  f^  *  0  that  the  growth  rate  of  instability  goes  through  two 
maxima  as  £  is  varied.  The  first  maximum  occurs  at  £  =“  1 ,  where  electro¬ 
static  effects  are  strong,  and  the  second  one  at  u>t  =“  1  where  the 

a 

instability  is  driven  by  the  magnetic  phase  lag.  For  a  highly  relativistic 
beam,  the  growth  rate  is  equal  at  the  two  maxima. 

In  Sec.  VI,  similar  models  and  analytic  techniques  are  applied  to 
the  theory  of  the  hose  instability  in  a  medium  of  arbitrary  conductivity, 
and  with  arbitrary  values  of  the  fractional  charge  and  current  neutraliza¬ 
tion.  The  instability  properties  are  similar  in  many  respects  to  those  of 
the  m  *  0  modes:  the  hose  instability  also  has  maximum  growth  rates  at 

£  88  1  and  wt  ~  1,  it  too  is  driven  strongly  unstable  by  return  current, 
d 

and  the  hose  growth  rates  are  comparable  to  those  of  the  sausage. 


4 


II.  EQUILIBRIUM  THEORY  AND  BASIC  ASSUMPTIONS 


The  equilibrium  configuration  consists  of  an  intense  electron  beam 
that  is  propagating  parallel  to  a  uniform  applied  magnetic  field 
through  a  background  plasma  with  a  conductivity  o (r) ,  and  is  located 
inside  a  grounded  cylindrical  conducting  wall  with  radius  R^.  Cylindrical 
polar  coordinates  (r,0,z)  are  used,  with  the  z-axis  along  the  axis  of 
symmetry.  In  equilibrium  (3/3t  =  0) ,  the  beam  is  assumed  to  be  azimuthally 
symmetric  (3/39  =  0),  infinitely  long,  and  axially  uniform  (3/3z  =  0).  The 
number  of  electrons  per  unit  axial  length  (N^)  is  defined  by 


<x>  0  2  0 

=  2ir  J  dr  r  n^(r)  =  TTR^n^r  =  0) 


where  n°(r)  is  the  equilibrium  density  profile,  and  R^  is  the  characteristic 
beam  radius.  In  the  present  analysis,  we  assume  that 


v 


me 


1, 


(1) 


2 

where  v  is  Budker's  parameter,  c  is  the  speed  of  light  in  vacuo,  Ybmc  is 
the  characteristic  electron  energy,  -e  and  m  are  the  electron  charge  and 
rest  mass,  respectively.  Equation  (1)  guarantees  that  the  beam  is  paraxial, 
i.e.,  transverse  electron  velocities  are  much  smaller  than  axial  velocities. 

In  this  paper,  we  investigate  the  equilibrium  and  stability 
properties  for  a  steady-state  (3/3t  =  0)  beam  distribution  of  the  form 


2irv.  m 
•b 


5  (H 


OJ 

b  9 


n  2 

-  yrac  )  6  (?z  -  Yhm  3bc) , 


f?(H,P.  ,PJ 
b  9  z 


(2) 


2  4  2  2  h 

where  the  total  energy,  H  =  (m  c  +  c  J3  )  -  e<|)0  ( r)  ,  the  canonical 

2 

angular  momentum,  P.  =  rp.  -  eBrtr  /2c,  and  the  axial  canonical  momentum, 
u  y  u 

Pz  =  Pz  -  (e/cjA^r) ,  are  the  three  single-particle  constants  of  the 

motion  in  the  equilibrium  fields-  Here,  p  =  (p  ,p.,p  )  is  the  mechanical 

%  r  o  z 

momentum,  n,  ,  w,  ,  and  y  are  constants,  y .  is  related  to 
b  b  b 

2  -% 

8,  by  y,  -  (1  -  6.  )  ,  A„  (r)  is  the  axial  component  of  vector  potential 

D  a  D  0 

for  the  equilibrium  azimuthal  self -magnetic  field,  and  <J>Q(r)  is  the 

equilibrium  electrostatic  potential.  As  a  consequence  of  Eqs.  (1)  and 

(2) ,  all  electrons  have  axial  velocities  V very  close  to  the  constant 

value  8.  c. 
b 

In  order  to  make  the  theoretical  analysis  tractable,  we  assume 
that  in  equilibrium  the  plasma  space  charge  provides  a  fractional 
neutralization  of  the  beam  charge  that  is  uniform  through  r  and  z,  i.e., 


0  0  0 
n.  (r)  -  n“(r)  =  f  n.  (r)  , 
i  e  e  b 


where  f  is  constant  with  0<  f  <1.  It  is  further  assumed  that  the 
e  —  e  — 

equilibrium  plasma  return  current  J^(r)  has  the  same  radial  profile  as  the 

P 

0 

beam  current  J,  (r) ,  i.e., 
b 

0  0 
J  (r)  *  -  f  J.  (r)  , 
p  mb 

where  f  is  a  positive  constant  with  0  £  f  £  1.  An  eigenvalue  treatment  of 
in  in 

the  instability  is  possible  only  if  f  and  f  are  taken  to  be  z  -  independent 

0  m 

(which  is  not  strictly  consistent  in  a  resistive  medium) .  Radial  profiles 
of  plasma  charge  and  current  which  differ  from  the  beam  profile  are  of 
interest,  but  are  not  studied  in  the  present  paper. 

Since  the  beam  is  paraxial,  it  is  straightforward  to  show  that  the 


term  H  -  w.  P.  in  Eq.  (2)  can  be  approximated  by" 

D  u 


(4) 


H  ■  Vb=  Ybmc 


px2 

2y7m  +  V0(r) 


2  2  2 

where  p  =  p  +  (p,  -Y,  muu  .  r)  is  the  transverse  momentum-scruared  in  a  ro- 
x  r  9  b  b 

tating  frame  with  the  angular  velocity  uu^,  and  the  effective  potential 
is  defined  by 


Vr)  =  iV(JJb 


'a,b)(JJb  -vr2' 


(5) 


In  Eq.  (5)  , 

,  2  A  2  1/2 

♦  >^a-y  -  a -*.)][  •  (.) 

where  uu  =  eBQ/Ybmc  is  the  beam  electron  cyclotron  frequency  and 
A  2  2 

=  4TTe  is  the  beam  electron  plasma  frequency-squared.  For 

radial  confinement  in  the  equilibrium  state,  it  is  required  from  Eq.  (6) 
that 


*2  2  2  A2 
01  3^(1  -  f  )  +  4)  .  /2  >  0)  .  (1  -  f  ) 
pb  b  m  cb  pb  e 


(7) 


which  assures  that  the  repulsive  space-charge  force  on  the  beam  electron  is 
weaker  than  the  magnetic  focussing  force. 

Making  use  of  Eqs.  (2)  and  (4),  we  find  the  beam  density  profile 


nb(r) 


nfa  ,  0  £  r  <  ^  , 

°  ,  ^  <  r  <  Rc  , 


(8) 


where  the  beam  radius  is  defined  by 

2  2  *  -1  +  -1  _  i 

■  2c  (y  -1)  (wb  -  uib)  )  •  T^e  equilibrium  solution  exists 
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for  rotational  frequency  satisfying  u>“  <  <  uu*.  Additional  equilibrium 

properties  associated  with  the  distribution  function  in  Eq.  (2)  are  discussed 
in  Refs.  3  and  5. 
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III.  LINEARIZED  VLASOV-MAXWELL  EQUATIONS 


In  this  section,  we  use  the  linearized  Vlasov-Maxwell  equations  to  ob¬ 
tain  an  integro-differential  eigenvalue  equation  for  the  azimuthally 
symmetric  (d/dQ  =  0)  modes  of  an  electron  beam.  We  adopt  a  normal  mode 
approach  in  which  all  perturbed  quantities  are  assumed  to  vary  with  z  and  t 
as 

A 

I  (x,t)  =2  (r)exp[  i(kz  -  uj  t  )]  , 


with  ImJ)  >  0.  Here  u)  is  the  oscillation  frequency  and  k  is  the  axial  wave- 

number.  For  paraxial  beams  with  all  electron  axial  velocities  approximately 

constant  at  v^  =  j3  c  [which  follows  from  (1)  and  (2)],  it  is  more  convenient 

to  use  T  and  z ,  rather  than  t  and  z ,  as  independent  variables",  where 

T  =  t  -  z/v  .  In  this  representation,  a  Fourier  decomposition  of  perturbed 
z 

quantities  is  expressed  as 

A 

2  (x,t)  =  I  (r)exp[-i(0JT  +  ftz/v  )]  ,  (9) 

<v  Z 


where  Q  =  0)  -  kv  =  uj  -  kq3  is  the  shifted  eigenfrequency  "seen"  by  a 

Z  .D 

beam  particle  (but  no  relativistic  charge  of  frame  is  implied) . 

The  subsequent  stability  analysis  is  restricted  to  a  long  wavelength 
and  low  frequency  perturbation  characterized  by 


IkRj2  «  1  , 


<< 


1  . 


c 


(10) 


If  we  forego  the  possibility  of  recovering  the  non-neutral  limit, 
a  o  ,  the  analysis  can  be  carried  through  with  the  weaker  assumptions 

|fi  Rj/cl  «  1  (11) 


9 


and 


!kVlrbl  ~  1  • 


(12) 


Equation  (11)  is  a  weak  form  of  "frozen  field"  assumption,  based  on  the 

fact  that  Cl  is  characterized  by  the  natural  oscillation  frequencies  of 

± 

particle  dynamics,  ui.  and  to  6.  ,  which  always  satisfy  (11)  in  a  paraxial 

°  pb  b 

beam.  Assumption  (12)  is  satisfied  either  for  long  wavelengths  or  in  a 
highly  relativistic  beam. 

The  perturbed  electric  and  magnetic  fields,  ^(x)  and  13  (x)  ,  can  be 
expressed  in  terms  of  perturbed  magnetic  and  electric  potentials  A(£)  and 
<t>(x)  : 


B  (x)  =  V  x  A(x)  , 

rw  /v  /v 


(13a) 


E(x)  =  i  -  A(x)  -  V<|>  (x) 


(13b) 


Introducing  the  Lorentz  gauge, 

V  .  A(x)  -  i  —  $  (x)  =0  , 

rs*  /v  »v  O  ^ 


(14) 


into  the  Maxwell  equations,  and  using  (11),  the  perturbed  potentials  in 
Eq.  (13)  satisfy 


(15) 


7x  ♦  =  -  4lTp 


(16) 
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where  p (x)  and  £(x)  are  the  perturbed  charge  and  current  densities, 
which  must  be  determined  self-cons istently ,  and  the  subscript  a  denotes 
transverse  components.  The  perturbed  current  density  contributed  by  the 
background  plasma  is  given  by 


J 

~P 


<&>  -  o<*)E(x)  *  a^i  ^A(x)  - 


(17) 


where  0 (r)  is  the  conductivity  of  the  background  plasma.  Defining  an 
effective  potential 


'l'(r)  =  A  (r)  -  <J>  (r) /6.  , 

z  b 


(18) 


and  using  (11) ,  the  axial  components  of  Eqs.  (15)  yield 


V2A(i  +  $/8b)  +  £ 

c  \  b 


(‘ 


s.n  ; 

c  bz 


(19) 


where  J  (x)  is  the  axial  component  of  the  beam  current  density. 

After  some  algebraic  manipulation  of  the  axial  component  of  Eq.  (13a) , 
Be  =7  x  A  .  we  can  show  that 

Z~Z  •"■'A  •"■'A 


cV  x  B  e  = 

~A  2~Z 


ikcV  *  +  i  |  0  -  + 

V  w 


(20) 


where  J,  is  the  transverse  component  of  the  perturbed  current  density, 

^  i3  a  unit  vector  in  the  z-direction,  and  use  has  been  made  of 

Eqs.  (15  -  18).  Making  use  of  Eqs.  (17)  and  (20),  it  is  also  straight¬ 


forward  to  show  that  for  axisymmetric  modes  (where  B  =  const) , 

z 


11 


1  3^_ 
r  3r 


(’ 


1 

r 


3_ 

3r 


ike  ( 3/3 r)^  -  4ttJ, 
_ _ _ br 

4iro-i(Q  -  kc/6.  Y?) 

d  b 


i  3_ 
+  r  3r 


(4tto(jj/c)A 

_ r _ 

4irc-i(n  -  kc/8 .  Y?) 

D  D 


(21) 


where  and  are  the  radial  components  of  the  perturbed  beam  current 
density  and  magnetic  potential,  respectively. 

A  JK 

From  the  continuity  equation  -iuip  +  £  .  j  =  0  ,  and  the  perturbed 
Poisson  equation,  it  is  evident  that  the  perturbed  electric  potential 

*  2 

<i>(r)  is  at  most  of  order  4ttR^V  .  vA/(4iro  -  io») .  Substituting  Eq.  (21)  into 
Eq.  (19)  and  using  (10)  we  therefore  obtain  the  eigenvalue  equation  for 
axisymmetric  modes. 


13  f  _ 4tl 

r  3r  )  r 

(  4TTO-1 


4  ita  -  itu 


i(Q  -  kc/BbYb 


fj  h  *(rl) 


4irao)i 


♦  (r) 


=  -4ir 


jbz(r)  +  1  3 

"brlr) 

c  6.  r  3r 

D 

4ira-i(ft  -  kc/B,Y^) 

D  O  . 

, (22a) 


where  the  perturbed  beam  current  density  J.  =  J,  e  +  J.  e  is  to  be 

''■©  bz'^z  br~r 

calculated  from  the  Vlasov  equation.  We  find  in  all  cases  that  J.  is  the 

DZ 

driving  term  in  the  regime  4iro>c,  but  that  must  be  retained  in  going 
to  the  small  a  limit* 

We  begin  our  analysis  of  the  Vlasov  equation  by  using  the  method  of 
characteristics  and  neglecting  initial  perturbations ,  to  write  the 
perturbed  beam  electron  distribution  in  the  form 


J°  &  exp(-inz/Bbc){E(x()  +  £  Lv’  x  B <&•)]}•  frr  ^  ' 

b 


12 


where  z  =  z'  -  z,  and  the  particle  trajectories  sc ’  (z1)  and  v'(z') 

***  rw 

satisfy  the  "initial"  condition  x'(z  =  0)  =  x,  and  v'  (z  ■  0)  =  v. 

#w  «v 

Within  the  context  of  Eq.  (10) ,  the  perturbed  beam  electron  distribution 
in  Eq.  (11)  can  be  further  simplified  as7 


ffa(r,p)  - 


_b  fb  ' 
x  3px  . 


iMr)  +  — 
v 

z 


0  ~  *  ~ 
p  dziiHr' )exp(-i\J>z/v  ) 


where  p^  is  the  transverse  momentum  in  a  frame  of  reference  rotating  at 

2  2  2 

frequency  i.e. ,  by  Px  =  Fr  +  (Pg  -  y^moi^r)  .  The  perturbed  beam 
current  density  can  then  be  expressed  as  a  velocity  moment  of  fb(r,£)  , 


4(r)  35  &J*  /  dpz£  {lx  ail 


i_rn_  !fb 

VzLPx  3px 


ymB.  c  3p 
o  z 


1  o  ^  , 

I  f*  ^ 

j  dzi)/»(r 

J  —00 


'Jexpf-ii/iz/v^) 


where  we  have  again  used  (10)  and  the  property  that  A^fr)  is  order  of 
4>(r)  or  less.  The  integro-differential  eigenvalue  equation  for  azimuthally 
symmetric  perturbations,  Eqs.  (22a,b) ,  constitutes  a  principal  result  of 
this  article,  and  can  be  used  to  investigate  stability  properties  quite 
rigorously  for  a  broad  range  of  system  parameters.  The  source  term  in 
Eq.  (22b)  contains  an  integral  of  the  unknown  eigenfunction  ip(r)  over  the 
equilibrium  particle  orbits,  which  makes  the  equation  rather  intractable, 
in  general.  However,  it  is  possible  to  obtain  analytic  solutions  to 
Eq.  (22)  for  some  ranges  of  physical  parameters.  In  the  remainder  of  this 
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paper,  we  investigate  the  stability  properties  of  the  analytic  solutions 


to  Eq.  (22) . 

In  order  to  carry  out  the  orbit  integration  in  Eq.  (22b) ,  it  is 

useful  to  introduce  Cartesian  coordinates  (x,y)  which  are  related  to  the 

polar  coordinates  (r,9)  by  x  =  r  cos  9,  y  =  r  sin  9,  and  to  note  that 

the  polar  momentum  variables  (p^,<J>)  in  the  rotating  frame  are  related  to 

p  and  p  by  p  +  Y.w.y  =  p  cos*,  and  p  -  y  mm  x  =  p  simji,  where 
x  y  x  bo  x  y  b  b  ± 

2  2  2 

p  =  p  +  (p.  -  y.mo).  r)  .  For  the  case  under  consideration,  where  the 
r  o  b  b 

unperturbed  net  current  is  uniform  out  to  r  =  R^,  the  perpendicular 
electron  trajectories  can  then  be  expressed  as^ 


x'(z)  = 


%  -  % 


-  {sin((j>  +  id.  z/v  )  -  sin($  +  o>z/v  )} 

Y,  m  0  2  0  2 


—  -H**  +  — *** 

+  r(ui^  -  o).)cos(9  +  id.  z/v  )  -  r(ui.  -  w.)cos  (9  +  u>,z/v  )>  ,  (23a) 

DO  02  OO  02 

y'  (z)  ■  -  < -  {cos(<J>  +  a),  z/v  )  -  cos(<f>  +  id,  z/v  )} 

+  -  1  y,  m  b  z  b  z 

“b  *  b 


+  r(w.  -  w^)  sin  (9  +  w,z/v  )  -  r(a>.  -  u,)sin(9  +  id  z/v  )f,  (23b) 

bb  bz  bb  d  z 


where  the  frequencies  <u,  are  defined  in  Eq.  (6).  In  the  next  section,  we 

o 

shall  develop  an  approximation  scheme  for  solving  the  coupled  equations 
(22),  based  on  exact  solutions  to  Eqs.  (22)  and  (23)  that  can  be  obtained 
in  the  limit  |  4iroR  id/c  |  «  1 . 


IV.  DISPERSION  RELATIONS  FOR  THE  AXISYMMETRIC  MODES 


We  shall  proceed  with  the  stability  analysis  of  the  azimuthally 
symmetric  modes  for  the  particular  configuration  illustrated  in  Fig.  1, 
wherein  the  plasma  conductivity  profile  is  specified  to  be  a  step 
function, 


a(  r) 


with  a2  so  small  that 


0<r<Rfa  , 

V  Vr<Rc 


(24) 


*  2  *  2 
°2  Rc  «  ai  *b 


(25) 


This  type  of  conductivity  profile  is  reasonable  for  a  beam  propagating 
in  and  ionizing  a  neutral  or  weakly  pre-ionized  gas.  The  dependence  of 
stability  properties  on  the  form  of  the  conductivity  profile  is  of  con¬ 
siderable  interest,  but  will  not  be  considered  in  full  generality  in  the 
present  paper. 

The  method  used  to  derive  approximate  dispersion  relations  is  as 
follows.  We  define  a  characteristic  magnetic  decay  time 

Td  =  ^1  2  "  (26) 

[Each  mode  has  a  different  characteristic  magnetic  decay  time  which 
differs  from  (26)  by  a  numerical  factor.  It  is  convenient  to  use  (26) , 
which  defines  the  dipole  decay  time  for  the  hose  mode,  as  a  reference 
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decay  time  throughout  the  paper.]  In  the  limit  tux  <<1,  solutions  of 

m 

Eqs.  (22)  can  be  found,  with  the  form 

n  . 

^(r)  =  £  a^r/I^)  3  ,  O^rs*^  ,  (27) 

j  =  o 

and  the  quadrature  of  Eq.  (22b)  can  be  performed  analytically  to  yield 
expressions  for  in  terms  of  ip .  On  the  assumption  that  these  expressions 
for  iMr)  are  reasonably  accurate  in  the  whole  range  0<ujt  <  1,  we  proceed  to 

A  A 

use  these  expressions  for  in  terms  of  t|i  in  the  full  Eq.  (22a)  ,  including 
the  oit^  term,  thus  reducing  Eqs.  (22)  to  an  ordinary  differential  equation 
for  0.  This  equation  could  be  solved  exactly  for  a  dispersion  relation 
S1(uj)  (in  terms  of  Bessel  functions),  but  in  the  spirit  of  the  preceding 
analysis  we  choose  rather  to  produce  a  variational  approximation  to  the 
dispersion  relation,  using  as  the  trial  eigenfunction.  One  or  more  modes 
are  found  for  each  radial  mode  number  n,  with  n  =  1  corresponding  to  the 
usual  sausage  mode  and  n  =  2  to  beam  hollowing  modes. 
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A.  n  =  1  Sausage  Mode 

As  a  first  example,  we  consider  the  n  =  1  eigenfunction-  In  the 
limit  of  small  magnetic  decay  time, 


3  “V<  1' 


(34a) 


the  eigenfunction 


-  ,  fo'^X1  '  0  sr  s  \ 

*V>  ,  lys  r  s  R. 


s  r  s  r 

c 


(35) 


is  a  self-consistent  solution  to  Eq.  (22)  ,  as  shown  in  the  previous 
literature.3  Since  we  perturb  about  this  solution,  the  stability  analysis 
in  this  section  is  restricted  to  the  regime 


-  (ox  <  1 
3  d  ~ 


(34b) 


The  perturbed  beam  current  density  can  be  calculated  in  closed  form  by 
substituting  Eq.  (35)  into  Eq.  (22b)  and  using  Eq.  (23) .  The  result  for 
Jbz,  the  axial  component  of  the  perturbed  current  density,  is 


Jbz(r) 


f!* 

4ir 


1  3_  ( 
r  3r  |3 


“L(r) 

_ 


a  -  («. 


-  v2 


(36) 


where  the  plasma  frequency  function  uj^(r)  is  defined  by 


A  2 
(  V 
)  °-  % 


0  s  r  < 


*b 


<  r  s  r 


C371 
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r 


and  the  fast  and  slow  rotational  frequencies  u),  for  beam  electrons  are 

b 

defined  in  Eq.  (6) .  Similarly,  the  radial  component  of  perturbed  current 
density  is  found  to  be 


3.  (r) 
br 


i 

4n 


Mpb(r)8bfl 


a2  - 


(ub  ' 


V 


3?*<r) 


(38) 


Equations  (36) -(38)  show  that  this  mode  consists  of  a  self-similar  expansion 
and  contraction  of  the  beam  slices ,  without  altering  the  flat-topped  beam 
profile  or  mixing  beam  electrons  at  different  axial  coordinates  r  within 
the  beam.  This  is  the  usual  sausage  mode. 

Substituting  Eq.  (36)  and  (38)  into  Eqs.  (22)  and  carrying  out  a  simple 
algebraic  manipulation,  we  obtain  the  eigenvalue  equation. 


-SIMM -  Lrc-iu, 

4ira-i(fl  -kc/8bYb)  L 

4iroaii  „ 

2  v(r)  =  0,  (39) 
c 

for  the  n  =  1  perturbation.  Equation  (35)  is  clearly  a  self  consistent 
solution  to  Eq.  (39)  in  the  limit  ut."*  0  (i.e.  the  last  term  can  be 

G 

neglected) ,  provided  that  ft(u>)  satisfies  the  dispersion  relation 


In  order  to  find  an  approximate  dispersion  relation  which  includes  the 
influence  of  the  magnetic  decay  time  without  going  through  a  Bessel 


_£b _ 


“O2 


H  *  A 


1  3 _ 

r  3r 


function  analysis  of  Eq.  (39),  we  multiply  Eq.  (39)  by  nMr)  and  integrate 


over  r  from  r  =  0  to  r  =  R  .  It  was  shown  in  Ref.  7  that  if  a  trial 

c 

function  is  substituted  in  the  integrals ,  this  procedure  gives  a 
dispersion  relation  (2(oj)  that  is  accurate  to  second  order  in  the  error 
in  i p  .  Thus  the  procedure  can  be  described  as  a  variational  approxi¬ 
mation,  although  it  does  not  give  a  lower  bound,  because  the  differential 
operator  is  non-Hermitian.  The  result  is  a  dispersion  relation 

A  *  W  ,  /  A  - 

(4ira1  -  iw)  (1  -  -  iurr^  +  - -  [4Tro  S  +  i  =  0  (40) 

a  ~  (a)b  "  V  '  Yb/ 

where  use  has  been  made  of  Eq.  (10) .  In  the  limit  of  a  nonneutral 
electron  beam  where  ^c^/ai  -*  0,  the  dispersion  relation  in  Eq.  (40)  can 
be  simplified  to 


2  +  -  2  A  2  2 
a  =  (wb  -  V  +V  ^b 


(41) 


which  is  identical  to  the  result  obtained  by  Uhm  and  Davidson.  The 


equilibrium  constraint,  Eq.  (7) ,  indicates  that  the  right-hand  side  of 
Eq.  (41)  is  always  positive.  Thus  the  n  =  1  perturbations  about  a  non¬ 
neutral  electron  beam  are  stable.  Analysis  of  Eq.  (40)  will  be  deferred 


B.  n  =  2  Hollowing  Mode 


As  a  second  example,  we  consider  the  n  =  2  eigenfunction. 


<Mr)  = 


ao(l-4r/R^  +  3r4/R^) ,  , 

0  ,  ^  ^  r  ^  Rc  , 


(42) 


which  is  found  to  be  a  self-consistent  solution  to  (22)  in  the  limit 


IF  "Td  ~  1 


(43) 


Substituting  Eq.  (42)  into  Eq.  (22b)  and  carrying  out  some  tedious 
algebraic  manipulations  that  make  use  of  Eq.  (10) ,  we  can  show  that  the 
radial  and  axial  components  of  the  perturbed  n  =  2  current  density  are 


Jbr(r) 


and 


i  “BhW8b°  U 

"  a2  -  <«r  -  <o2  t9r 


<Mr) 


'b  ~b' 


.  VVV^i  %(2.3  % 

a2  -  4(^  -  «-)2  R^\  R2 


(44) 


Jbz(r) 


c6 

b  1  3_ 
4ir  r  3r 


<d2.  (r) 

-PL 


.2  ,  +  -2  3r 

n  -  (ub  ~  V 


fr  <Hr) 


cBb  72a0(Mb :  %){% :  v 


4tt  r  2  .  +  •".2-7  r  J2  .  .  +  “ .  2  - 

Lft  -  ((i>b  -  ajfe)  J  [ft  -  4(u>b  -  0)b)  J 
(■  ■  ■  ?) 


<5  (r  - 


V 


Pb  ^ 


UDb(r) 

*4JS- 


(45) 


where  is  the  rotational  frequency  of  the  electron  beam.  We  note  that 


2 

J,  (r)  is  of  the  form  A  +  A, r  within  the  beam,  where  A  and  A,  are 
bz  o  1  o  1 

constants.  Thus  this  mode  involves  non-self-similar  axisymmetric 

distortions  (hollowing  or  peaking  on  axis)  of  the  beam  profile. 

As  in  the  previous  section,  an  approximate  eigenvalue  equation 

for  wt.tO  can  be  determined  by  substituting  Eqs.  (44)  and  (45)  into 
a 

Eq.  (22) .  This  gives 


1  3_  (  / _ 

r  3r  (Wo  -  i 


4ira-iw 


i((2  -  kc/BbYb) 


>)  'h  *<rl| 


2  +  Fb(r) )  r  fr  ^(r)[+  <Hr) 


=  72a, 


(wb  '  V  (“b  "  “fc 


°  ft2  -  4(ol 


rb,r)£2 

U  - 3  4V 

< 

vJ 

(46) 


where  the  function  T,  (r)  is  defined  by 

D 


rb(r)  *  - 


«L(r) 

pb 


2  2 

4-rro6.  +  iw/y. 

b  b 


«.  +  -  2  2 
n  -  (u),  -  to,  )  4irc-i(S)  -  kc/B.Y.)  • 

D  D  DO 


(47) 


As  in  the  previous  section,  we  multiply  Eq.  (46)  by  nj>(r)  and  integrate 

over  r  from  r  =  0  to  r  =  R  ,  to  obtain  a  variational  expression  for  the 

c 

dispersion  relation  in  terms  of  the  unknown  exact  i|i(r)  ,  and  then  use 
Eq.  (42)  as  the  trial  function.  This  procedure  gives 


(4ira1 


+ 


a 


<4,"’i6b 


Is  l“v) 


+  iu/Yb)^ 


18(u,b  -  a>b)  (a>b  -  U‘) 

2  +  -  2 
fl  -  4  (a),  -  a),  ) 
b  b 


(48) 


where  use  has  been  made  of  Eq.  (10).  The  dispersion  relation  in  Eq.  (48) 
can  be  used  to  investigate  stability  properties  of  the  n  =  2  perturbations 


for  a  broad  range  of  physical  parameters,  and  detailed  investigations 


of  Eq.  (48)  will  be  presented  in  Sec,  V. 

As  a  particular  case,  we  consider  Eq.  (48)  in  the  limit  of  a 
nonneutral  electron  beam  characterized  by  4ttc^/w  -*  0.  In  this  case, 
the  dispersion  relation  in  Eq.  (48)  can  be  expressed  as 


r  _2  4*  —  2-i  r  *2  +  —  2-1 

[ft  -  (id  -  id  )  J  Lft  -  4  (id,  -  id,  )  J 
b  b  b  b 


-f  {[ft2 


4  (id,  -  Id.  )  ]  +  18  (id  -  id,  )  ((d,  -  Id.  )  } 

D  D  D  D  D  D 


which  is  identical  to  the  results  obtained  by  Uhm  and  Davidson.3 

2 

Equation  (49)  is  a  simple  quadratic  equation  for  ,  and  the  necessary 
and  sufficient  condition  for  instability  can  be  expressed  as 


(“b  "  a>b) 

,  +  2 
(«db  -  <db) 


m’  ■  ■] 


Note  that  when  Eq.  (50)  is  satisfied,  the  perturbations  are  purely 

growing,  i.e.  ,  (2  =  Reft  =  0,  and  that  the  rotational  frequency  co,  can 

r  b 

have  a  large  influence  on  stability  behavior.  Equation  (49)  is  also 

16 

similar  in  form  to  the  result  obtained  by  Gluckstem  for  transverse 

instabilities  of  proton  beams  in  the  quadrupole  magnetic  fields,  which 

12  13 

are  particularly  important  in  the  heavy  ion  fusion  experiment.  ' 


r 


C.  n  =  0  Axial  Bunching  Mode 

The  n  =  0  mode,  which  is  of  particular  interest  for 
unneutralized  beams,  is  quite  different  in  nature  from  the  n  ^  1  modes. 
We  find  that 

(1,  0  S  r  <  R  , 

*(r>  -  ao 

Un(r/Rc)/An(Rb/Rc)  ,  <  r  £  Rc  (28) 

is  a  solution  to  Eq.  (22)  in  the  limit 


4uit.  <<  1  . 
d 


(29a) 


and 


*b  K< 


R  . 
c 


(29b) 


Substituting  Eq.  (28)  in  Eq.  (22b)  and  using  (29b)  and  the  identity 
Pr  =  p^cos  (<p-8)  ,  we  obtain  the  axial  component  of  the  perturbed  beam 
current  density, 


(r) 


2"  2  3  2 

-  e  nbaQw  /(YbmcQ  )  , 

0  , 


0  *  r  <  Rb  , 
*b  <  r  *  Rc  ' 


(30) 


and  the  radial  component. 


Jbr (r)  =  °  ’ 


(31) 


Clearly,  the  n  =  0  mode  involves  purely  axial  flow,  leading  to  bunching. 
Because  a  beam  becomes  axially  rigid  in  the  highly  relativistic  limit. 


is  seen  to  fall  off  rapidly  as  increases.  The  dispersion  relation 


that  determines  the  eigenfrequency  a>  or  52  is  obtained  as  in  the 
previous  sections,  by  substituting  Eqs.  (28),  (30),  and  (31)  into 
Eq.  (22)  ,  multiplying  Eq.  (22)  by  rtp(r)  and  integrating  over  r  from 
r  =  0  to  r  =  R  .  This  gives  the  dispersion  relation 


id  +  4iria2 

[(kc/B^)  -  id  -  4ifia2  ]ta  (R^/R^) 


+  4itor 


d 


where  v  is  Budker's  parameter  defined  in  Eq. 
only  if 


(1), 


(32) 


a  result  which  is  valid 


4<DTd  SI,  (29c) 

and  (29b)  is  also  satisfied. 

The  dispersion  relation  in  Eq.  (32)  can  be  used  to  investigate 
stability  properties  of  the  azimuthally  symmetric  surface  perturbations 
about  an  electron  beam  in  a  background  plasma  with  arbitrary  values  of 
density  and  plasma  conductivity.  In  particular,  for  a  nonneutral  electron 
beam  in  an  environment  where  4ttct/<d  “♦  0,  Eq.  (32)  simplifies  to 


n2  2v  n2  2 

a  -  —  (k  c 
Yb 


-  <d  )£n(Rc/Rfa) 


(33) 
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which  is  identical  to  the  result  obtained  by  Briggs  for  the  space-charge 

wave  mode  where  oo  *■  kcS^-  Note  from  Eq.  (33)  that  the  space-charge  wave 

2  2  2 

is  a  stable  mode  because  k  c  >  u  in  this  mode.  We  also  note  the 
-3  2 

dependence  of  52  ,  which  occurs  because  this  mode  (unlike  the  nil 
modes)  involves  purely  axial  flow.  A  detailed  investigation  of  Eq.  (32) 


will  be  carried  out  in  Sec.  V. 


V. 


STABILITY  PROPERTIES  OF  AZIMUTHALLY  SYMMETRIC  PERTURBATIONS 


We  now  investigate  stability  properties  predicted  by  Eqs.  (32) , 
(40) ,  and  (48)  for  a  %  0.  Depending  upon  the  way  in  which  a  stability 
problem  is  phrased,  either  u>  or  0  may  be  regarded  as  the  independent 
variable,  and  taken  to  be  real;  a  complete  understanding  requires  an 
analysis  in  the  complex  plane  for  both  oi  and  ft.  We  shall  concentrate 
here  on  the  solution  of  the  dispersion  relations  in  Eqs.  (32)  ,  (40)  , 
and  (48) ,  for  given  real  values  of  the  oscillation  frequency  w;  then  ft 
determines  the  growth  of  the  wave  on  a  particular  beam  segment  as  the 
beam  propagates  downstream. 


1 

l 


f 

f 
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A.  Stability  Analysis  in  the  High  Conductivity  Regime 

To  make  the  theoretical  analysis  tractable,  we  assume  that 
the  plasma  is  collisional  to  the  extent  that  it  is  characterized  by  a 
real,  scalar  conductivity.  Moreover,  we  consider  in  this  section  the 
case  in  which  the  perturbed  beam  space  charge  field  is  completely 
neutralized  by  the  plasma,  which  requires  a  high  conductivity  plasma, 

4itc  «  to  .  (51) 


Within  the  context  of  Eq.  (51) ,  the  dispersion  relations  in  Eqs.  (32) , 
(40) ,  and  (48)  can  be  simplified  considerably. 

(a)  When  (51)  is  satisfied,  the  dispersion  relation  in  Eq.  (32) 
for  n  =  0  reduces  to 


2v  to  . . 

-  -  -  4l»t  - 

Tb  fl 


ta<W 


(52) 


which  clearly  indicates  instability.  If  in  addition  the  magnetic  decay 
time  is  small,  |  4uiTdJin  (R^/R^)  |  <<  1,  Eq.  (52)  simplifies  further  to 


n2  =  -  ~  to2)ln(Rc/Rb)  ,  (53) 

which  yields  a  purely  growing  mode.  Note  from  Eq.  (53)  that  the  growth 

3/2 

rate  Q.  is  inversely  proportional  to  y  ,  which  indicates  that  this 
l  b 

unstable  mode  is  important  only  for  a  mildly  relativistic  electron  beam 
satisfying  y  <  5;  the  rapid  fall-off  of  (2.  with  y.  is  due  to  the  purely 

D  ^  ID 

axial  flow  that  occurs  in  this  mode.  It  is  instructive  to  compare 
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Eq.  (53)  with  Eq.  (33),  in  which  the  terms  proportional  to  k  c  and 
2 

oi  are  due  to  the  perturbed  beam  space  charge  and  magnetic 
fields,  respectively.  In  a  nonneutral  electron  beam,  the  stabilizing 
influence  of  the  perturbed  space  charge  field  always  prevails  over  the 
destabilizing  influence  of  the  perturbed  magnetic  field;  thus  the  beam 
supports  only  stable  space  charge  waves.  However,  in  a  high  conductivity 
plasma  characterized  by  Eq.  (51) ,  the  stabilizing  influence  of  the 
perturbed  space  charge  field  vanishes;  the  mode  is  then  unstable,  as 
shown  ir  Eq.  (53) . 

(b)  The  dispersion  relation  in  Eq.  (40)  for  n  =  1  simplifies 


-2  "2  2 
•  V  *  “pbsb 
,  +  -.2 
{ub "  V 


3  laiTd  ' 


when  4irc  »  o.  The  instability  mechanism  can  be  easily  identified  in  a 

self-pinched  electron  beam  where  the  applied  magnetic  field  vanishes  and 

the  fractional  charge  neutralization  becomes  unity  (oi  ,  =  0,  f  =1). 

cb  e 

In  this  particular  case,  Eq.  (54)  can  be  expressed  as 


2  "2  2 

8  -  </.  fT(l  -  2  f  ) 

po  o  m 


2  "2  2 
n  -  210.6.(1  -  f  ) 

pb  o  m 
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Equation  (55)  clearly  indicates  that  the  instability  of  the  n  *  1 


sausage  perturbation  is  driven  by  the  plasma  return  current  (f  )  and  the 

m 

magnetic  phase  lag  (got^)  .  Even  for  a  very  small  magnetic  decay  time 

satisfying  <<  1,  the  system  supports  a  purely  growing  mode  whenever 

the  partial  current  neutralization  f  satisfies 

m 


(56) 


f  >0.5; 
m 


the  growth  rate  of  this  mode  is 


Q. 

1 


VBb(2V1) 


>5 


We  note  further  that  Eq.  (55)  reduces  to 


z2  -  n  -  32  [i  +  i  |  (a  Vc)2^  '  (57) 

2 

for  the  case  uuT  ,  <<  1.  In  the  limit  of  ri  =  2j3.  ,  Eq.  (57)  is  remarkably 
a  b 

g 

similar  to  the  result  obtained  by  Weinberg  within  the  framework  of  the 
coaxial  circular  electron  beam  model. 

The  influence  of  the  magnetic  decay  time  on  the  stability  behavior 
will  be  investigated  numerically  in  Sec.  VB. 

(c)  The  dispersion  relation  in  Eq.  (48)  for  n  =  2  simplifies  to 


18uj  2J,2  to*  -  ujJ  to,.  -  uO 


pb  b 


[Q  -(au; 


— .  2 *ip  ^2 

UUb)  J[u  - 


4(*b  '%)2] 


A  O2 

(U  ,3. 

pb  b 


o2  ,  +  — ,  2 

Q  -  to.  -  oj,  ) 
b  b 


(58) 


-  -1  *  i?  im  d 


for  high  conductivity  plasmas.  We  note  from  Eq.  (58)  that  for  n  *  2  as 

well  the  magnetic  phase  lag  (ujt^  )  and  plasma  return  current  (fm)  cause 

the  instability.  In  order  to  illustrate  the  influence  of  the  plasma 

return  current  on  the  stability  behavior,  we  consider  Eq.  (58)  for  a 

self-pinched  beam  (uj  ,  =  0,  f  *  1)  and  for  a  small  magnetic  decay  time 

cb  e 
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( out  <<  1) .  In  this  particular  case,  the  dispersion  relation  in  Eq. 
a 

(58)  can  be  further  simplified  as 

Z4  -  [10(1  -  f  )  -  1]Z2  +  [16(1  -  f  )+  1J  (1  -  f  )  *  0  ,  (59) 

m  m  m 

A 

where  the  normalized  eigenfrequency  Z  is  defined  by  Z  =  Q/u  3  ,  and  we 

pb  Jo 

have  also  set  ou.  =  0,  i.e.  we  consider  the  case  of  a  non-rotating  pinched 
b 

2 

beam.  Equation  (59)  is  a  simple  quadratic  equation  for  Z  ,  and  the 
necessary  and  sufficient  condition  for  instability  can  be  expressed  as 


f  >  (2  -  3  '  /2)/3  =  0.38  .  (60) 

m 

In  contrast  to  the  case  of  n  =  1,  we  note  from  Eq.  (59)  that  the  n  =>  2 
perturbations  are  not  always  purely  growing  when  Eq.  (60)  is  satisfied, 
i.e.,  Z^fO.  Unstable  n  «  2  perturbations  can  reduce  the  beam  density 
at  axis,  leading  to  a  hollow  beam  profile.  Detailed  numerical  investi¬ 


gations  of  Eq.  (58)  are  carried  out  in  the  next  section. 


B.  Numerical  Analysis  for  Stability  Properties 

In  this  section,  we  summarize  the  results  of  numerical  studies 
of  the  dispersion  relations  in  Eqs.  (32),  (40),  and  (48)  for  a  broad 
range  of  the  beam  electron  energy  y.  ,  the  normalized  oscillation  frequency 

D 

jur^/c,  the  plasma  conductivity  parameter 


C  =  4TTa]/u>  , 


and  the  focussing  parameter 


+  -  2  A  2 

n  -  (^b  -  »b)  /*pb 


(62a) 


The  analysis  in  this  section  is  restricted  to  a  self-pinched  electron  beam 

with  oj  *  0.  [The  stabilizing  influence  of  the  applied  magnetic  field 
cb 

(ju  .  )  on  a  related  instability  (resistive  hose)  is  discussed  in  Ref.  7.] 
cb 

Making  use  of  Eq.  (6) ,  we  can  simplify  the  parameter  T|  in  Eq.  (61)  to 

t)  =  2[^(1  -  f  5  -  (1  -  f  )]  ,  (62b) 

D  in  e 

for  self-pinched  beams.  We  note  that  the  value  of  the  parameter  T]  is 

2 

limited  to  the  range  0  £  r|  £  2j3 

Defining  the  normalized  eigenfrequency 

Z  =  C2/uu  ,  (63) 

pb 

the  dispersion  relations  in  Eqs.  (32),  (40),  and  (48)  can  be  expressed 


2  *"(yy  /*% 
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w,/^)C3^  +  i  Aft 


(*%\  (g2/gl^b  + 

V  C  )  [(a  /aL)C  -  ill  -  i 


£  (uiRjj/c)  jin(Rc/Rb)] 


for  the  n  =  0  perturbation. 


Z  -  r\  + 


2  2 
*3b  + 


(£  -  i)[l  -  iC  tui^/c)  /6] 


=  0 


(65) 


for  the  n  =  1  perturbation,  and 

2  2 

2  2  t  1  '■/Yw 

(Z  -  T))  (Z  -  4T|)  +  (Z  +  -  n)  - 5 - - - 5 -  =  0  , 

<£  -  i)[l  -  iC  «JUR/c)Z/30] 

for  the  n  =  2  perturbation.  In  obtaining  Eq.  (64)  ,  use  has  been  made  of 
Eq.  (10).  Moreover,  we  have  assumed  uu.  =  0  in  Eq.  (66),  which  is  con- 

D 

sistent  for  self-pinched  beams.  Equations  (64) ,  (65) ,  and  (66)  clearly 
show  that  only  the  n  =  0  mode  is  influenced  by  the  location  of  the  con¬ 
ducting  wall. 

Typical  numerical  results  of  Eqs.  (64)  -  (66)  are  summarized  in 

2 

Fig.  2  for  Yb  =  1.5,  (wR^/c)  =  0.01  and  the  parameter  rj  «  1.11.  Note 

from  Eq.  (62)  that  rj  =  1.11  is  the  largest  possible  value  of  T|  for 

Yb  *  1.5,  corresponding  to  a  space  charge  neutralized  electron  beam 

(f  =1)  with  no  plasma  return  current  (f  =  0) .  We  note  the  following 
e  m 

features  of  Fig.  2:  (i)  The  n  =  0  and  n  =  1  unstable  modes  have  Z  <  0, 

r 

i.e.  propagate  backward  in  the  beam,  while  one  of  the  two  n  =  2  unstable 
modes  [denoted  n  =  2(+)]  has  Z^  >  0  and  the  other  [n  =  2(-)]  has  Z^  <  0. 
The  n  *  2(-)  mode  is  always  the  faster  growing  of  the  two.  (ii)  Each 
mode,  except  n  =  0,  has  two  maxima  of  the  growth  rate  £1^,  the  first 
occurring  for  all  modes  at  £  »  1,  the  transition  point  from  electrostati¬ 
cally  to  magnetically  dominated  forces,  and  the  second  occurring  at 

'  a»  efc/DR^)  2,  i.e.  JUT  as  3/4,  (n  =  1)  , 


(66) 


(67a) 


Note  that  when  f  «  1,  the  second  instability  peak  is  driven  by  the  phase 
m 

lag  T  between  the  beam  oscillations  and  the  magnetic  restoring  force,  and 

it  is  appropriate  to  scale  uu  to  T^.  (iii)  The  n  =  1  perturbation  is  the 

most  unstable  azimuthally  symmetric  mode.  (iv)  For  a  space-charge  neutralized 

beam  (f  =1)  with  f  <<  1,  the  absolute  value  of  Doppler-shifted  real  fre- 
e  m 

1  /0  1  /? 

quency  can  be  approximated  by  |z  |  *•  r\  for  n  =  1,  and  |z  |  Tj  and 

i  i  1/2 

\Z^\  2r\  for  n  =  2  over  the  entire  unstable  range  of  the  conductivity 

parameter  £  [Fig-  2(b)].  Since  this  dependence  is  so  simple,  Z^_  will  not  be 

shown  in  the  subsequent  figures. 

Figure  3  is  similar  to  Fig.  2a,  except  that  a  higher  frequency 
2  2  2 

to  R^  /c  =  0.1)  is  considered.  We  note  that  Eqs.  (67)  for  the  location 
of  the  magnetic-dominated  instability  peak  are  well  satisfied;  the  location 
of  the  peak  thus  moves  to  the  left  on  a  £  scale.  The  first  instability  peak 
remains  at  £  1.  The  growth  rate  of  the  n  =  0  mode  is  seen  to  increase 

significantly  with  toR^/c) ,  but  the  n  =  1  and  n  =  2  modes  depend  weakly  on 
this  parameter. 

In  Fig.  4  we  consider  a  highly  relativistic  (yfa  =  10)  case,  again  with 

f  *  1,  f '  *  0  (which  gives  t)  =  1.98  in  this  case) .  The  n  =  0  mode  is  not 

shown  since  its  growth  rate  is  very  small  for  large  y,  .  The  most  striking 

b 

change  from  Figs.  2  and  3  is  that  the  two  peak  growth  rates  of  each  mode 

have  become  equal,  and  in  fact  the  Z(£)  curves  are  symmetric  (on  a 

2 

logarithmic  scale)  about  a  point  £  ~  toR^/c)  .  This  feature  is  easily  seen 
to  be  a  consequence  of  Eqs.  (64)  -  (66).  When  y^  >>  1,  each  of  these  equa¬ 
tions  takes  the  form 
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F  (Z) 


const 


(68) 


(i  -  Aj/C>  (i  -  a2/C) 

where  F(Z)  is  a  function  of  Z  that  is  different  for  each  of  the  two  inodes, 
and  A^  and  A2  are  constants  which  take  different  values  for  each  of  the 
two  modes.  It  follows  that  Z(Q)  takes  the  same  value  at  two  points  and 
C2  =  (A^/A2£ 1) ,  i.e.  the  Z(£)  curve  is  symmetric  (on  a  log  plot)  about 
C  =  (A^)172. 

The  dependence  of  stability  properties  on  the  parameter  T]  is  illus¬ 
trated  in  Fig.  5(a)  for  Yb  =  1.5,  n  =  1,  in  Fig.  5(b)  for  yfa  =  1.5, 

n  =  2,  and  in  Fig.  5(c)  for  *  10,  n  =  1 ,  where  the  normalized  growth 

2 

rate  Z^  is  plotted  versus  C,  for  (uiP^/c)  =  0.1  and  several  values  of  the 
parameter  T).  It  is  evident  that  as  the  parameter  T)  [defined  in  Eq.  (62)] 
decreases  to  zero  (e.g.  weakly  pinched  beam-near ly  complete  current  and 
space  charge  neutralization) ,  the  growth  rate  increases  rapidly,  and  the 
strong  growth  occurs  over  a  broad  range  of  u; ,  rather  than  being  confined 
to  two  sharp  peaks.  As  discussed  in  Sec.  V.A  for  the  high  conductivity 
regime,  this  remarkable  behavior  is  due  to  repulsion  of  the  beam  current 
by  the  plasma  return  current  (fm  >  0)  flowing  in  the  highly  conducting 
plasma  channel.  We  also  note  from  Fig.  5(b)  that  the  n  =  2  (+)  perturbation 
is  stable  for  rj  =  0:  the  dispersion  relation,  Eq.  (66) ,  reduces  to 

+  i/Yb _ 

(C  -  i)C  1  -  K  toRj/c)  2/30] 

which  gives  only  one  unstable  mode,  n  *  2  (-) . 
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VI.  HOSE  INSTABILITY 


Although  this  paper  has  been  concerned  primarily  with  the  axisymmetric 
(m  =  0)  instabilities,  essentially  the  same  analysis  can  be  used  to  treat 
higher  m  modes  over  the  full  range  of  values  of  the  background  conductivity 
O.  provided  that  J(r,T)  is  r-independent  over  the  range  0  <  r  <  R^.  In 
particular,  it  is  of  interest  to  compare  the  growth  characteristics  of  the 
m  =  1  hose  instability  to  those  of  the  axisymmetric  modes  we  have  studied. 
Our  analysis  of  the  hose  instability  will  be  based  on  the  same  equilibrium 
{Sec.  II)  except  that  we  now  take  a  to  be  a  constant,  which  can  have  any 
value,  over  the  full  range  of  r  and  2.  We  need  not  make  the  long  wavelength 
assumption  (10);  the  weaker  assumptions  (11)  and  (12)  suffice. 

For  the  hose  mode,  all  perturbed  quantities  are  taken  to  vary  as 


4(x,t)  =  <J>  (r)  exp  i(kz  +3  «-  uut) 


(69) 


=  <|>(r)  exp  l-Kuut  +  Ciz/v  -  0)]  , 


so  that 


?2  =  1_  d_  a _ 

x  r  dr  r  dr  2  dr  r  dr  r 


The  analysis  of  Maxwell's  equations  follows  that  of  Sec.  II  through  Eq.  (19), 
at  which  point  the  constancy  of  a  is  used  in  (18)  and  (19)  to  show  that 
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The  analysis  of  the  Vlasov  equation  is  similar  to  that  of  Sec.  II,  leading 
in  this  case  to  the  expression 
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+  (0  -  ii),  )  j  dz  if  (r"*)  exp  [  i (6  9  -  rfz) ] 


Equations  (70)  define  the  m  =  1  eigenvalue  problem  as  Eqs.  (22)  did  for 

m  *  0. 

In  the  case  of  the  hose  mode,  in  the  limit 
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Eqs.  (70)  support  a  solution 
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essentially  a  distortion  free  snake  like  displacement  of  the  beam  and  the 

associated  electromagnetic  fields  (within  the  constraints  imposed  by  the 

conducting  boundary  conditions  at  r  =  R  ) .  Using  the  orbit  equations  (23) , 

c 

Eq.  (70b)  can  then  be  solved  for  the  perturbed  current,  with  the  results 
Wr)  "  ■  Tn  - l^L6(r-B),  (73a) 
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(73c) 


^b0(r>  =  1  Jbr(r)  +  ^  Jbz{r)- 

b 


For  finite  values  of  out  we  follow  the  same  approximation  scheme 
used  in  Sec.  IV  to  derive  the  dispersion  relation.  Equation  (73)  is  used 
in  Eq.  (70a)  to  derive  the  dispersion  relation. 
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We  note  that  in  the  high  conductivity  limit  4TTO  >>  |uu  |  ,  Eq.  (74)  reduces 
to  the  eigenvalue  equation  previously  obtained  by  the  authors,7 
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As  in  Sec.  IV,  we  proceed  by  multiplying  both  sides  of  Eq.  (74)  by  r  y  and 
integrating  over  r  from  r  «*  0  to  r  =  R^,  thereby  obtaining  a  variational 
expression  for  the  dispersion  relation  Q(uu), 

-  J’  ^  ^4  to  -  iai  + 


*'  dr  /._  .  ^b>“rt.2V> 

2  (fl-aj. ' +)  (n-i u~) , 

D  D 


5?  <*♦>] 


~4  n  i  am 


(4  to  -  iuu  )  j  c  drr  £ 


A  2 


(76) 


uu 


4rra  -i  I  w  - 


kc  _  pb 

Vb2 "  3b2 


A  <5  ■)  A  4 

*jA 


(f2-0)b  )  (£l-u»  ) 


36 


We  then  use  Eq.  (72)  as  the  trial  function  (r)  in  (72),  finally  arriving 


at  an  approximate  dispersion  relation 


4TTJ 
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Equation  (77)  can  be  used  to  investigate  hose  stability  properties  over  a 

broad  range  of  system  parameters,  including  conductivity  J,  fractional 

charge  and  current  neutralization  f  and  f  ,  applied  maanetic  field  strength 

e  m 

(vJUcb)  ,  electron  energy  (yfa)  ,  and  proximity  to  the  conducting  guide  (R^/R^)  . 
In  the  limit  of  4na  <<  |uu  j  of  low  background  conductivity,  Eq.  (77)  reduces 
to  the  result  obtained  by  Uhm  and  Davidson,^ 


gf  = 


*  2 
JJ 

.  .  Eb- 


2y.  2(fi-ou  ")  (fl-w. +) 

b  b  b 


(78) 


which  leads  only  to  stable  oscillations.  On  the  other  hand,  in  the  high 


conductivity  regime  4tt<j  >>  |cu  |  ,  Eq.  (77)  reduces  to  our  previous  result, 
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For  comparison  with  the  results  of  Sec.  V  for  the  axisymmetric  in¬ 
stabilities,  we  present  a  number  of  numerical  solutions  of  Eq.  (77)  for 
Re  ft(ju)  and  ImQ(dU),  specifying  j)b  =  0  and  =  0-5,  but  sur¬ 

veying  a  broad  range  of  beam  electron  energy  yb/  normalized  frequency 

oaR.  /c,  normalized  conductivity  C,  ,  and  focussing  parameter  r|.  These  re- 
o 

suits  have  been  plotted  as  the  dashed  curves  in  Figs.  2-4,  along  with 
the  results  for  the  various  axisymmetric  modes  (solid  curves) .  The  hose 
dispersion  relation  is  seen  to  be  remarkably  similar  to  those  of  the 
axisymmetric  modes.  In  Figs.  2-4  (cases  with  no  return  current) ,  the 
hose  growth  rate  n.  is  seen  to  have  two  sharp  peaks,  the  first  at  Q  «  1 
and  the  second  at 

2 

'  «  3  (c/U)  Rfa)  ,  i.e.  jJTd  «  0.4  .  (80) 

Both  the  location  of  the  peak  growth  rate  given  by  Eq.  (80)  and  the  value 

of  the  peak  growth  rate  depend  to  some  extent  on  the  value  of  ;  the 

corresponding  results  for  the  axisymmetric  modes  are  independent  of  R  . 

c 

The  peak  hose  growth  rates  plotted  in  Figs.  2-4  for  ■  §  lie  between 

those  of  the  n  =  0  and  n  =  1  axisymmetric  modes;  for  /Rc  -*  0,  however, 
the  hose  growth  rates  would  be  slightly  greater  than  those  of  the  axisym¬ 
metric  modes.  In  the  limit  of  large  y  ,  Eq.  (77)  can  be  put  in  the  form 
(68),  so  that  the  two-peaked  structure  of  Z(£)  becomes  symmetric  on  a  log 
C  plot,  as  seen  in  Fig.  4  for  all  modes.  As  the  return  current  increases, 
i.e.  I)  decreases  to  zero,  the  growth  spectrum  of  the  hose  mode  broadens  and 
becomes  singly  peaked,  and  the  maximum  growth  rate  increases  significantly. 
Like  the  sausage  and  hollowing  modes ,  the  hose  is  driven  strongly  unstable 
by  the  mutual  repulsion  between  the  beam  current  and  return  current,  when 
the  latter  has  the  same  equilibrium  radial  profile  as  the  former. 
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VII.  CONCLUSIONS 


In  this  paper,  we  have  investigated  the  stability  properties  of  both 
azimuthally  symmetric  O/d0  =  0,  or  m  =  0)  and  hose  (m  =  1))  perturbations 
in  an  intense  particle  beam  propagating  parallel  to  a  uniform  applied 

A 

magnetic  field  B(^ez  through  a  background  plasma.  The  analysis  was  carried 

out  within  the  framework  of  the  Vlasov-Maxwell  equations,  assuming  long 

wavelength,  low  frequency  perturbations.  The  analysis  was  simplified  by 

considering  only  step- function,  time- independent  conductivity  profiles, 

by  requiring  that  the  plasma  charge  and  current  density  have  the  same 

radial  profile  as  the  beam,  and  by  considering  a  beam  with  a  "loss  cone" 

distribution  function,  in  which  all  of  the  beam  electrons  have  the  same 

value  of  axial  canonical  momentum  and  the  same  value  of  energy  in  a  frame 

of  reference  rotating  with  angular  velocity  uu.  .  This  distribution  leads 

b 

to  a  flat-topped  beam  radial  profile;  the  unperturbed  beam  electron  orbits 
are  then  simply  harmonic  at  a  single  frequency. 

Equilibrium  properties  were  calculated  in  Sec.  II.  The  formal  sta¬ 
bility  analysis  for  azimuthally  symmetric  perturbations  was  carried  out  in 
Sec.  Ill  and  an  integro-dif ferential  eigenvalue  equation  was  obtained, 
which  included  beam  electron  thermal  effects  and  applied  generally  to  any 
radial  mode  number.  In  Sec.  IV,  dispersion  relations  for  the  axisymmetric 
modes  with  radial  mode  numbers  n  *  0,  1,  and  2  were  obtained  analytically 
from  the  integro-dif ferential  eigenvalue  equation  (22) ,  for  a  moderate  value 
of  the  magnetic  decay  time  such  that  uut^  ^  1.  Also  discussed  in  Sec.  IV  were 
the  stability  properties  of  a  nonneutral  electron  beam  characterized  by 

A 

4 TT »  0.  A  parallel  treatment  of  the  hose  mode  was  carried  through  in 
Sec .  VI . 
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Analytic  and  numerical  investigations  of  the  n  =  0,  1,  and  2  axi- 


symmetric  dispersion  relations  were  carried  out  in  Sec.  V,  and  the  hose 

dispersion  relation  in  Sec.  VI,  for  a  collision-dominated  background 

plasma.  It  was  shown  that  the  n  *  1  and  n  =  2  axisymmetric  modes  and 

the  hose  mode  are  all  driven  unstable  by  both  the  plasma  return  current 

(f^)  and  the  magnetic  decay  time  (jut^)  •  Even  for  a  very  small  magnetic 

decay  time,  the  system  can  be  unstable  whenever  the  fractional  current 

neutralization  satisfies  f  >0.5  for  n  =  1  and  f  >0.38  for  n  *  2. 

m  m 

We  found  that  the  n  =  1  perturbation  is  the  most  unstable  axisymmetric 

mode,  with  maximum  growth  rate  *  0.205,  for  a  relativistic  beam  with 

f  *  0f  slightly  less  than  that  of  the  resistive  hose  instability  when 

R^/R^  =  0,  but  greater  than  the  hose  growth  rate  when  ^  g.  Finally, 

it  has  been  found  in  the  case  of  zero  plasma  return  current  (f  =0)  that 

in 

the  growth  rate  of  instability  has  two  local  maxima  corresponding  to  the 

A 

electrostatic  regime  (^tror^  •  uu )  and  to  the  regime  where  the  instability 

is  driven  by  the  magnetic  decay  time  (UUT^  ~  1) . 

The  resistive  hose  instability,  in  the  high  conductivity  regime,  has 
6  7 

been  studied  '  for  rounded  beam  profiles  which  introduce  dispersive 
anharmonic  effects.  It  is  of  considerable  interest  to  extend  the  analysis 
of  the  axisymmetric  modes  in  this  way,  as  well  as  to  study  the  effects  of 
more  general  conductivity  and  plasma  current  geometries.  It  is  believed 
that  the  hose  and  axisymmetric  modes  respond  differently  to  these 
generalizations  of  the  equilibrium.  Such  effects  will  be  considered  in 
future  work. 
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Fig.  1  -  Plasma  conductivity  profile  [Eq.  (24)]. 
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